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Transverse energy and charged particle pseudorapidity densities at midrapidity and their ratio, 
dEr/dr) \ m id /dN c h/dn\ m id, are evaluated in a statistical model with longitudinal and transverse 
flows for the wide range of colliders, from AGS to RHIC at ^snn = 200 GeV. Evaluations are 
done at freeze-out parameters established independently from fits to observed particle yields and 
Pt spectra. Decays of hadron resonances are treated thoroughly and are included in derivations of 
dEr I drj\ m ici and dN c h/ dn\ m id- The predictions of the model agree well with the experimental data. 
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However, some (explicable) overestimation of the ratio has been observed. 
C3 . PACS numbers: 25.75.-q, 25.75.Dw, 24.10.Pa, 24.10.Jv 



I. INTRODUCTION 



(N 

In this paper, the idea of an independent test of the applicability of a statistical model for the description of the 
soft part of particle production in a heavy-ion collision postulated in [l| , is developed for the much more realistic case 
■ of a hadron gas and its expansion. So far, the statistical model has been applied successfully in description of particle 
yield ratios and px spectra measured in heavy-ion collisions B,ISI^i np,p ,Ml^ lj^lj^lj^ll^lj^lT5| (there are also 
computation packages for thermal studies available from the Web [la, ll3)- Now, the freeze-out parameters obtained 
from those analyses will be used to evaluate global observables: the transverse energy density dEx/dr/, the charged 
particle multiplicity density dN c h/drj and their ratio. The advantage of such an approach is based on the fact that 
transverse energy measurements are independent of hadron spectroscopy (in particular, no particle identification is 
Q-i! necessary), therefore they could be used as an additional test of the self-consistency of a statistical model. The same 
holds true for the charged particle multiplicity, which actually is the charged hadron multiplicity, according to the 
| experimental definition given in |l8j . 

The experimentally measured transverse energy is defined as 

> : 

^ ■ L 

E T = ^E i -sm6 i , (1) 

d ; i=i 

where 0i is the polar angle, Ei denotes Ei — {uin means the nucleon mass) for baryons and the total energy Ei 
for all other particles, and the sum is taken over all L emitted particles 0|. Additionally, in the case of RHIC at 
y/SNN = 200 GeV, Ei + rriN is taken instead of Ei for antibarvons I2CM . 

The statistical model with single freeze-out is used (for details see Fyjfand references therein). The model reproduces 
very well ratios and px spectra of particles measured at RHIC |J llOi ll lj . The main assumption of the model is the 
simultaneous occurrence of chemical and thermal freeze-outs, which is important if px spectra are considered (this 
enables to neglect the possible elastic interactions after the chemical freeze-out). Since in the present paper the 
integrated quantities over px are dealt with, the above-mentioned assumption should not be so important for final 
results. 

The actually detected (stable) particles have two sources: (a) a thermal gas and (b) secondaries produced by decays 
and sequential decays of primordial resonances. All stable hadrons and confirmed resonances up to a mass of 2 GeV 
from the Particle Data Tables [2l| are constituents of the gas. The distributions of particles from source (a) are given 
by a Bose-Einstein or a Fermi-Dirac distribution at the freeze-out. The distributions of secondaries (source (&)) can 
be obtained from the elementary kinematics of a many-body decay or from the superposition of two or more such 
decays (for details see the Appendix and 0). In the following, all possible (2-, 3- and 4-body) decays with branching 
ratios not less than 1% are considered. Also almost all possible sequential decays are taken into account, namely: 
2o2, 2o2o2, 2o2o2o2, 2o3, 2o4, 3o2, 3o3, 2o2o3, 2o3o2, 3o2o2, 2o3o3, where 2, 3 and 4 mean the 2-, 3- and 
4-body decay respectively, and a cascade proceeds from the right to the left (as in the usual mathematical definition 
of the superposition of functions). It should be stressed that all contributions from weak decays are included. The 
contribution to the transverse energy from the omitted cascades has been estimated at 0.2%. But since the most of 
these cascades ends with two photons, they do not contribute to the charged particle multiplicity at all. 
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II. THE BASIS OF THE FREEZE-OUT MODEL 

The following are the foundations of the model. A noninteracting gas of stable hadrons and resonances at chemical 
and thermal equilibrium is created at the Central Rapidity Region (CRR) of a collision. The gas cools and expands, 
and after reaching the freeze-out point it ceases. The conditions for the freeze-out are expressed by values of two 
independent thermal parameters: the temperature T and the baryon number chemical potential ps- The strangeness 
chemical potential ps is determined from the requirement that the overall strangeness of the gas equals zero. 

A freeze-out hypersurface is defined by the condition 



t 2 — r% — Ty — r| = const , (2) 

which means that the freeze-out takes place at a fixed moment of the invariant time r. Additionally, it is assumed 
that the four-velocity of an element of the freeze-out hypersurface is proportional to its coordinate, 



,M — _ - - 



T T 

Then the following parameterization of the hypersurface is chosen: 



t = Tcoshan coshaj_, r x = r sinh a_i_ cos 4>, r v = r sinh ax sin cf>, r z = r sinh ay cosh a_i_, (4) 
where a\\ is the rapidity of the element, an = tanh _1 (r 2 /i), and a±_ determines the transverse radius 



r x + r y — T sinh a^. (5) 

To keep the transverse size finite, r is restricted by the condition r < p max - In this way one has two additional 
parameters of the model, r and p m ax, connected with the geometry of the freeze-out hypersurface. 
Also the transverse velocity, v p , can be obtained 



r y \ 2 t sinh a± tanh aj 



t J \ t J t cosha|| cosha|| 



(6) 



which is the value of the transverse velocity /3j_ from the central slice after boosting it in the longitudinal direction. 
The transverse velocity can be expressed as a function of the transverse radius 

(3 x (r) = tanh ox = , = = • (7) 
V t + r 

Since it is an increasing function of r, the maximum value of (3±_ called the maximum transverse-flow parameter (or 
the surface velocity), is given by 

pmax _ Pmax _ Pmax/ T ^ 

V /r2 + Pmax V 1 + (Pmax/r) 2 ' 

so it depends only on the ratio Pmax/i~. 

III. TRANSVERSE ENERGY AND CHARGED PARTICLE DENSITIES 

According to the general description founded in and developed in 0, ^| for the case with decays taken into 
account, the invariant distribution of the measured particles of species i has the form 



d 2 p T dy 



p^da^frip-u), (9) 
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where da^ is the normal vector on a freeze-out hypersurface, p ■ u = p^u^ , is the four-velocity of a fluid element 
and fi is the final momentum distribution of the particle in question. The final distribution means here that h is the 
sum of primordial and simple and sequential decay contributions to the particle distribution (for details see 13]). For 
the hypersurface and expansion described in sect. [HI eq. 10 takes the following form: 

r 3 / dan / sinho;_L(i(sinhQ!j_) / d£ p- u fi(p-u) , (10) 



d 2 p T dy 



— ao 

where 



p ■ u = wlt cosh ay coshau — Pt cos£sinha_L ■ (11) 

Note that the distribution expressed by eqs. I|10[) and (|11|) is explicitly boost invariant (in fact, it is constant with 
respect to rapidity). 

The rapidity density of particle species i is given by 



dNi [ , 2 

whereas the corresponding pseudorapidity density reads 



dN t 
d 2 p T dy 



(12) 



dNi 



dq 



dy dNi 



drj d 2 pT dy 



d Pt 77- 



P dNi 



Ei d?p T dy 



Analogously, the transverse energy pseudorapidity density for the same species can be written as 



(13) 



= / d 2 PT Ei 



dr/ 

For the quantities at midrapidity one has 



p T dy dNi 



p drj d 2 pT dy 



= I d 2 p T pt 



E % dNi 



Ei d 2 p T dy 



(14) 



dNj 
dr] 



mid 



d pt 



dNi \Jp T + v c. m . s mi 



d 2 pT dy 



rriT 



(15) 



dE 



T,i 



dr) 



d 2 p T dy 



-, i = nucleon 



raid 



J d 2 p T p T d /J^ z dy , i + nucleon 



(16) 



where v c , m . s is the velocity of the center of mass of two colliding nuclei with respect to the laboratory frame (only for 
RHIC v c , m .s = 0). For RHIC at ^snn = 200 GeV the case i ^ nucleon in eq. (|To)) is replaced by 



<JN, 



dE T< 
drj 



m T +m N ■ _ 



antinucleon 



I d2 PT PT 3^7^ "^T 

J d 2 pT Pt d^p^dy ' ^ ^ nuc ^ eon ^ antinucleon 
Now, the overall charged particle and transverse energy densities can be obtained 



mid 



(17) 



dN cl 
dr) 

dE T 
dr) 



mid 



mid 



E 



i£B 



dNi 
dr] 



E 



dE Tii 
drj 



mid 



mid 



(18) 



(19) 



where A and B (B C A) denote sets of species of finally detected particles. In the view of the definition given in 
|l8j and the detailed description of the experimental setup and the analysis procedure from [2^, the set of charged 
particles B can consist of stable hadrons only, B = {tt + , n~ , K + , K~, p, p}, whereas A also includes photons, 
K° L , n and n |H|. 
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TABLE I: Comparison of the statistical model estimates of the rapidity densities of charged particles with the experimental 
values for the 5% most central collisions at -^/sjvjv = 130 GeV at RHIC l24| . 

Particles dN c h/dy\ y= o 



Theory Experiment 

7T+ + 7T" 548.6 546 ± 54.5 

K + + K~ 84.6 87.2 ± 12.0 

p + p 55.8 48.8 ±6.2 



IV. RESULTS 



To check the self-consistency of the described model, rapidity densities of pions, kaons, protons and antiprotons 
have been calculated for the 5% most central Au-Au collisions at y / s]v r /v = 130 GeV at RHIC. The supplied values of 
the thermal and geometric parameters are in this case T = 165 MeV, ps — 41 MeV, r = 8.2 fm and p m ax = 6.9 fm 
0, . The geometric parameters were obtained from the fit to the pr spectra of the above-mentioned particles [24| . 
The integrated yields over pr are also given in [24j , so the comparison with the predictions of eq. l|T^|) can be done 
easily. The results are presented in TABLEQ] Note that the very good agreement has been found. 

The presentation of the main results of the paper needs a few comments concerning the AGS case. In the all cited 
papers the same method of establishing the thermal parameters T and ps is applied. The method is based on the 
best fit of calculated particle density ratios to the relative particle abundance data. But the different model of the 
freeze-out was applied for the description of pr spectra measured at AGS |2j,|5(. In that model (for details see |25J), 
the freeze-out happens instantaneously in the r direction, i.e. at a constant value of t (not at a constant value of 
t as here). The shape of a hypersurfacc in the longitudinal direction is not determined explicitly, but due to the 
factorization of the transverse mass spectrum it can affect only the normalization. The parameters connected with 
the expansion are the surface velocity f3™ ax and p max - The transverse velocity profile has the following form 



0x(r)=(3T ax ( — ) , (20) 

\ Pmax / 

with the choice a = 1. Therefore, the implementation of values of (3™ ax obtained within that model into the presented 
one is entirely ad hoc , ne vertheless it works surprisingly well. Of course, one directly could apply the description of the 
transverse flow from |25j to calculate the transverse energy and charged particle densities, but it is much more tempting 
and elegant to work within one model. Additionally, there is one technical problem connected with the treatment of 
resonance decays. Here, the very convenient form of the invariant distribution, eq. @, has been derived because the 
normal vector is proportional to the four- velocity, da^ oc u^. This is not the case for the hypersurface chosen in 1251 . 
so the calculation of resonance decay contributions would be much more complex (for the exact formulae, see [l3jL 
Note also that in the view of |2t| . particle distributions depend very weakly on the exact form of the velocity profile 
(i.e. for considered a — 0.5, 1, and 2) in the model described in |25j. It can be checked that profile (JJJ lies in between 
two profiles of the form l|20|l with a = 0.5, and 1. 

To put values of /3™ aa: from 0,0 into formulae of sect. IIIII one should invert eq. (jSJ to obtain 



n arnax 

Pmax HJ_ \ 

t ~ yi - (f3™ ax y ' 

It should be recalled here, that the value of r itself is not necessary to calculate the transverse energy per charged 
particle, since this parameter cancels in the ratio. 

The final results of numerical estimates of dEj- / drj\ mi( i and dN ^ / drj\ m i,i together with the corresponding experi- 
mental data are listed in TABLElTll To make predictions for the AGS case it has been assumed that the maximal 
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TABLE II: Values of dEr / dr/\ m id and dN c h / dri\ m id calculated in the framework of the statistical model with expansion. In 
the first column thermal and geometric parameters are listed for the corresponding collisions. In the third and last column 
experimental data for the most central collisions are given. 



Collision case 


dE T /dr]\ 


rmd [GeV] 




/drj\ m id 




Theory 


Experiment 


Theory 


Experiment 


Au-Au at RHIC at ^sjviv = 200 GeV: 
T = 165.6 MeV, b = 28.5 MeV 
Pma* = 7.15 fm, r = 7.86 fm (f3^ ax = 0.67) [12] 

Au-Au at RHIC at ^s NN = 130 GeV: 
T = 165 MeV, p B = 41 MeV 
p max = 6.9 fm, r = 8.2 fm {f3^ ax = 0.64) [13J 

rb-Fb at bra: 
T = 164 MeV, fi B = 234 MeV 
p max = 6.45 fm, r = 5.74 fm (p?** = 0.75) [14, 15] 


585 
507 
447 


597 ± 34 [20] 
503 ± 25 [19] 
363 ± 91 [27] 


589 
555 
476 


699 ± 46 [20] 
579 ± 29 b 
[30] 

622 ± 41 [23] 
568 ± 47 b 
[24] 

464+f 3 [2JJ 


Au-Au at AGS: 
T = 130 MeV, fiB = 540 MeV 

prnax = Q g^ ^ = g 53 fm ^ J 


224 


» 200 [28] 


271 


« 270 [29] 


Si-Pb at AGS: 
T = 120 MeV, fiB = 540 MeV 

prnax = q = 5 Q2 fm gj 


57 


« 62 [29] 


91 


w 115 - 120 
[29] 



"For the modified definition of By, i.e. Ej + mjy is taken instead of Ej for antibaryons, see eq. iTi . 
6 For the charged particle multiplicity expressed as the sum of integrated charged hadron yields. 



transverse size p m ax equals the average of radii of two colliding nuclei and the nucleus radius has been expressed as 
Ra = tqAs , ro = 1.12 fm. Generally, the overall agreement is good. For RHIC the 11% — 16% underestimation of 
the charged particle density has been received (5% — 10% with respect to the lowest allowed values). But this result 
simply reflects the existing inconsistency in measurements of the charged particle multiplicity at RHIC. Namely, the 
sum of integrated charged hadron yields (see TABLEP), after converting to dN c h/dr] |2fJ, is substantially less then the 
directly measured dN c h/ 'drj\ m id- This is shown explicitly in the last column of TABLEllTI For RHIC at y / sjvw = 130 
GeV, the sum is 8.7% smaller then the total charged particle multiplicity. For RHIC at ^/snn — 200 GeV it is even 
worse, the sum is about 17% below the total dN c h/ drj\ mi d- But both values of the sum agree very well with the model 
predictions. Since the geometric parameters were established from the fits to the same px spectra, the agreement had 
to be obtained. Also for AGS the results agree qualitatively well with the data, in spite of the roughness of the method 
applied for this case. The overall error of evaluations of transverse energy and charged particle densities is about 0.5% 
and has two origins: (a) omission of the most complex cascades; (6) simplifications in numerical procedures for more 
involved cascades. The velocity of the center of mass of two colliding nuclei, v c . m . s , equals: for RHIC, 0.994 for SPS 
Pb-Pb collisions at 158 ■ A GeV, 0.918 for AGS Au-Au collisions at 11 • A GeV and 0.678 for AGS Si-Pb collisions at 
14.6 • A GeV. 

Values of the ratio dEx / 'drj\ m id / 'dN c h / 'dr]\ m id can be also given. They are collected in TABLE lTTTI together with the 
corresponding data. Generally, the overall overestimation of the order of 15% has been obtained. In the RHIC case this 
is the result of the underestimation of dN c h/ 'dr]\ m i4, which has been explained earlier. But when in the denominator 
of the experimental ratio, dN c h / drj\ mi( i from the summing up of integrated hadron yields is put, the theoretical 
predictions agree very well with the data. Note that the similar inconsistency in charged particle measurements 
could have also been the origin of the discrepancy between model and experimental values of dN c h/drj\ m id, seen in the 
AGS Si-Pb case. For SPS, the result agrees with the experimental value within errors. The overall error of model 
evaluations of the ratio is less than 1%. These results have been also depicted together with the data in fig-G] One 
can see that the relative positions of theoretical points agree very well with the data, they are shifted up only and 
this is the effect of the overestimation discussed earlier. 

It should be stressed that dEx I 'drj\ m id / 'dN ~ c h / 'dr]\ m id depends substantially on the value of f3™ ax . It can be seen 
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TABLE III: Values of the ratio dET/dri\ m id/dN c h/dri\mid calculated in the framework of the statistical model with expansion. 
In the last column experimental data for the most central collisions are given. 

Collision case dE T /dri\ mid /dN ch /dr]\ m i d [GeV] 



Theory Experiment 



Au-Au at RHIC at ^/s]7n = 200 GeV 0.99 " 0.87 ± 0.06 [20j 

1.03 ±0.08 b 

Au-Au at RHIC at v /s^7 = 130 GeV 0.91 0.81 ± 0.06 [19] 



0.89 ±0.09 
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Pb-Pb at SPS 0.94 0.78 ± 0.21 [2JJ 

Au-Au at AGS 0.83 0.72 ± 0.08 [29] 

Si-Pb at AGS 0.63 0.52-0.54 [29J 

"For the modified definition of Et, i.e. Ej + mjv is taken instead of Ej for antibaryons, see eq. IT1 , 

'Author calculations with the use of experimental values given in TABLE|TT1 and the denominator expressed as the sum of integrated 
charged hadron yields. 



TABLE IV: Dependence of the transverse energy per charged particle on the maximum transverse-flow parameter fl™ ax for 
Au-Au collisions at AGS. 

PT X dEr/dri\ mid /dN ch /dr)\ mid [GeV] 

0.4 0.61 
0.57 0.72 

0.675 0.83 
0.8 1.07 



from TABLE ITvI where the transverse energy per charged particle estimates have been listed for a few values of (3™ ax 
for Au-Au collisions at AGS. Having compared with the experimental data (see TABLE llTTll . one can notice that this 
model yields the value of j3™ ax which is slightly lower than the value obtained within the model described in pij . 



V. COMPARISON WITH A STATIC CASE 



It would be very interesting to check how expansion influences the transverse energy per charged particle. The 
expansion produces additional energy so this process should increase the energy of a particle emitted from a thermal 
source. The preliminary analysis of a static case was done in Q] . But to compare with the present results up-to-date 
calculations should be performed. In a gas with only 40 species (including antiparticles) was examined and feeding 
charged particles from weak decays of neutral resonances was excluded. Thus, to extract the expansion contribution 
to dEr I dr}\ m id/ dN c h/ dr]\ m id one has to apply the general scheme of sect . 11111 again, but with the proper replacement 
of the invariant particle distribution. 

For a static gas {static in the c.m.s, of course), the invariant distribution of the measured particles of species i has 
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FIG. 1: Values of the transverse energy per charged particle at midrapidity for the most central collisions. Black dots denote 
evaluations of the ratio in the framework of the present model (the second column of TABLE lllll . Also data points for AGS 
[29t (a circle for Au-Au and a vertical bar for Si-Pb), SPS (triangle), RHIC at ^/sTTiv = 130 GeV 193 (square) and RHIC 
at ^/sjvjv = 200 GeV |2(» (star) are depicted. For RHIC, points with the sum of integrated charged hadron yields substituted 
for the denominator are also depicted (crosses). 



the form 0] 



dN,. 



d 2 pT dy 



V* E* ME* 



(22) 



where E* is the c.m.s energy of the zth particle and V* denotes the c.m.s volume of the gas at the freeze-out. Thus, 
at midrapidity, one has: 



dN t 



d?pr dy 



V* m T fi(m T ) 



(23) 



Now the general formulae of eqs. H15|) - (|19fl can be applied, but with dNi/d 2 prdy given by eq. (|23|l instead of eq. I |10| l. 
The results of numerical evaluations of dEx/drj | mid/dN^/dr] \ m id for the static gas are collected in TABLElVI In this 
case only two (thermal) parameters are needed and they are the same as in TABLElTTl Having compared with 
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TABLE V: Values of the ratio cLEt /dri\ m id/dN c h/drj\ m id calculated for the static gas. In the last column experimental data 
for the most central collisions are given. 



Collision case 



dE T /dr]\ mid /dN ch /dri\ mid [GeV] 



Theory 



Experiment 



Au-Au at RHIC at ^/JWn = 200 GeV 
Au-Au at RHIC at v /s^T = 130 GeV 

Pb-Pb at SPS 

Au-Au at AGS 

Si-Pb at AGS 



0.89 " 
0.82 
0.71 
0.62 
0.53 



0.87 ±0.06 [20J 
0.81 ±0.06 [19] 
0.78 ±0.21 [2JJ 
0.72 ±0.08 [29J 
0.52-0.54 f29] 



"For the modified definition of Et, i.e. Ej + mjv is taken instead of Ej for antibaryons, see eq. IT! . 



TABLE IHIl one can see that expansion is responsible for the following increases of the transverse energy per charged 
particle: 11% for RHIC, 32% for SPS, 34% for AGS Au-Au collisions and 19% for AGS Si-Pb collisions. This can 
be explained reasonably. The transverse energy per charged particle has two contributions: the first thermal and the 
second originated from expansion. The first is governed mainly by the temperature and the second by the maximum 
transverse- flow parameter (3™ ax . For a given temperature, the increase of j3™ ax should cause the weighting of the 
expansion contribution. But for a constant value of /3™ aa; , the strengthening of this contribution can be maintained by 
the lowering of the temperature. This is why for almost the same [3™ ax (see the first column of TABLE lTT|) the relative 
growth of the transverse energy per charged particle, after switching the expansion on, is much greater for AGS Au-Au 
collisions than for RHIC ones. On the other hand, for comparable temperatures, the expansion contributes to the 
transverse energy per charged particle much stronger for SPS than for RHIC (the former has substantially greater 



VI. CONCLUSIONS 



The expanding thermal hadron gas model has been used to reproduce transverse energy and charged particle 
multiplicity pseudorapidity densities and their ratio measured at AGS, SPS and RHIC. The importance of the present 
analysis originates from the fact that the transverse energy and the charged particle multiplicity are independent 
observables, so they can be used as new tools to verify the consistency of predictions of a statistical model for all 
colliders simultaneously. The predictions have been made at the previous estimates of thermal and geometric freeze- 
out parameters obtained from analyses of measured particle ratios and pr spectra at AGS SPS 0,0] and 
RHIC pjl The overall good agreement, not only of the ratio but also absolute values of dEr/dn \ m id and 
dN c h/ 1 dr\\ m idi with the data has been achieved. And the observed discrepancies can be explained reasonably. This 
strongly supports the idea that the thermal expanding source is responsible for the soft part of the particle production 
in heavy-ion collisions. Moreover, the description of various observables is consistent within one statistical model. 

In fact, there are additional arguments which make the above statement even more valuable. In principle, one 
could think at first glance that this analysis is nothing more like a kind of an internal consistency check of various 
measurements. And such a check could be done even in an model-independent way simply by integrating spectra of 
stable particles (the first time with the expression for transverse energy to obtain dET/drj \ m id and the second time 
without, to receive dN c h / drj^id) and then adding them all. But there are two reasons that this can not be done 
without any external input. First, transverse momentum spectra are measured in limited ranges, so very important 
\ow~pt regions are not covered by the data. For instance at RHIC, the first point for pions is at pr = 0.25 GeV/c, for 
kaons at pr = 0.45 GeV/c and for protons and antiprotons at pr = 0.65 GeV/c |24|,|3(j. There are also upper limits, 
but contributions from ranges above them are suppressed strongly in comparison with the low-px regions. Therefore, 
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to obtain integrated yields some extrapolations below and above the measured ranges are used. Usually two functions 
are used for each species and the contributed value is the average of their integrals. In fact these extrapolations are 
only analytical fits without any physical reasoning, but, for instance, contributions from regions covered by them 
account for 30% of the yield for pions, 40% for kaons and 25% for protons and antiprotons for RHIC at ^/snn = 130 
GeV 24]. On the other hand, a calorimeter acts very effectively for these species in the low-p-r range, namely pions 
with pt < 0.35 GeV /c, kaons with pt < 0.64 GeV/c and protons and antiprotons with pr < 0.94 GeV/c deposit all 
their kinetic energy |l9| . Since the very accurate predictions for the transverse energy density at midrapidity have 
been obtained (see TABLE lTT|) . the present analysis can be understood as an undirect proof that in these unmeasurable 
Pt regions spectra are also explicable by means of the thermal source with flow and decays. 

Second, it is impossible to check the consistency of the transverse energy data because not all stable hadron spectra 
are measured at midrapidity for each collision case. This mainly concerns neutrons and K^. The lacking contribution 
from hadron decay photons could be approximated to some extent with the use of ir° and r\ spectra, but they are 
also limited in ranges. And again, the very good agreement of model estimates of the transverse energy density at 
midrapidity with the data can be interpreted as the strong argument that the production of neutral stable particles 
can be described in terms of the expanding thermal source with superimposed decays. 

And last but not least, in opposite to the transverse energy, there is some inconsistency (of the order of 10%) of 
the independent measurements of charged particle multiplicities with the corresponding sums of integrated charged 
particle yields at RHIC (see sect. IIVf) . However, only for the case of ^/JWn — 200 GeV the substantial gap (of the 
order of 6% with respect to the direct measurement) between error bars of these two differently obtained values of 
dN c h I 'dr)\ m id exists. For the case of ^s^m — 130 GeV the error bars overlap almost one half of each other. But since 
the data at ^/s^n = 200 GeV are still preliminary, it is difficult to judge whether this inconsistency has the physical 
or experimental (an additional systematic error?) reason. 

The role of expansion is substantial. It produces about 10% — 30% of the transverse energy per charged particle. 
But, as can be seen from TABLElVl the expansion is not necessary to explain the experimental data for dEx/drj \ m id 
/dN c h/drj\ m id- The results suggest that the most of the transverse energy per charged particle is produced by 
the thermal movement. For sure, the expansion is necessary to explain the absolute values of dEx/dr) \ m id and 
dN c h/di~i\ m id- To obtain these one needs a volume of a place of "action" and the most adequate way to do it is to 
parameterize the evolution of the system in space and time, that is to put the expansion in. 

As the last comment, it should be stressed that the results of the present paper have been obtained within the model 
where the chemical freeze-out happens simultaneously with the thermal one. However, so far the most extensively 
studied scenario is that where the thermal freeze-out occurs later then the chemical freeze-out (for a review, see 
[3l| and references therein). This problem has not been addressed here. But one should notice that the distinction 
between these two freeze-outs means the introduction of the next parameter (the fifth here) into the model. Of course, 
an extra parameter in a phcnomcnological model always causes (or at least should cause) better agreement with the 
data. At the present level of investigations both spectra (refs. 0, 0> 0) an d global observables dEx/d-q \ m id 
and dN c h/dr/\ m id (this analysis) are predicted accurately with the assumption of the one freeze-out. However, more 
detailed studies should be performed to check whether the transverse energy per charged particle measurement could 
help somehow in distinction or not between these two freeze-outs and this will be the subject of further investigations. 
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APPENDIX 



The derivation of the momentum distribution of a product of a two-body decay M — ► m\ + mi can be found in 
|32fl or 0. For an n-body decay M — ► mi + mi -I- ... + m n , the momentum distribution of the product (labeled 1) 
can be written (following the method presented in 9jj) as: 

f[ n \\q],M, mi ,m 2 ,..., mn ) = B 2 * M + * 1 fd 3 kf M (\k\) 

2si + 1 iVW(M;mi,m2, —,m n ) J 
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// n j3 -> \ n n 

Vi=l 1 / i=l i=l 

x^Cigft-g), (A.l) 

where 

%ft = ft + { (7* - 1)^ + § } * , (A.2) 



7fc = lf ' E M = ^M* + k\ E i = y /m^+f i , (A.3) 

and sm (si) is the spin of the resonance (the product), B is the branching ratio, /m(|^|) denotes the momentum 
distribution of the decaying resonance and N^- 71 ' is the corresponding phase-space integral: 

N^(M; mi ,m 2 ,...,m n )= / M] <5(M - ^^ 3 )(^ft) . (A.4) 

\i=l 4 / i=l i=l 

The invariant amplitude for the decay, Ai, is assumed to be a constant here, so \Ai\ 2 cancels during normalization. 

With the use of the well-known technique of splitting up the phase-space integral into a convolution integral over 
two phase-space integrals (here, the first responsible for the 2-body decay and the second representing the (n— l)-body 
decay) [3^ |. the following recursive formulae for the n-body decay can be derived: 

f[ n \\q\,M,m l7 m 2 , ...,m„) = B — — - 



2s x + 1 Ar(n)(M;mi,m 2 ,...,m„) qEi(q) 

M—mi 

dm mN^ 71-1 ' (to; TO2, to„) 



where 



m.2 + ...+m n 

k-y- (g; Af,mi ,m) 

dkkf M (k), (A.5) 

k- (g; Af,mi ,m) 
M—mi 

A r(n) (M;TOi,TO 2 , ...,to„) = y dm m A r(rl_1) (m;m 2 , ...,to„) p(M;toi,to) , (A. 6) 



M 

fc±(g;M,TOi,TO 2 ) = — |p(M;TOi,TO 2 )£; 1 (g)± gJ B(M;mi,TO 2 )| , (A.7) 



for mi 7^ 0, whereas 



k + (q; M, toi, m 2 ) = +oo, 

|i(M 2 -ml) 2 -A*Vl 



k_(q; M,m 1 ,m 2 



(A. 



(M 2 -m 2 )g 



for mi = 0, and 



M (2) y/[M^ - (TO! + gg - (TO! -l^p] 

p(M;m 1 ,m 2 ) = —N ( >(M;m 1 ,m 2 ) = — 



(A.9) 
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E(M; mi ,m 2 ) = M2 -^ + m ' . (A.10) 
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